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Micromechanism  Based  Modeling  of  Structural. 
Life  in  Metal  Matrix  Composites 


David  H.  Allen  and  Dimitris  C.  Lagoudas 

Center  for  Mechanics  of  Composites 
Texas  A&M  University 
College  Station,  TX  77843-3141 

ABSTRACT 

AFOSR  Grant  F49620-94-1-0341 

This  report  details  the  most  recent  accomplishments  that  have  lead  towards  the  fulfillment  of  the  grant 
objectives.  These  achievements  include:  1)  life  prediction  of  continuous  fiber  metal  matrix  composites;  2) 
the  influence  of  heat  treatment  on  the  mechanical  properties  and  damage  development  in  a  SiC/Ti-15-3 
MMC;  3)  the  experimental  characterization  of  oxidation  on  fracture  surfaces  and  crack  growth  behavior;  4) 
modeling  the  effects  of  oxidation  on  the  crack  growth  resistance  of  metals;  and  5)  the  modeling  of 
oxidation  fronts  in  metals. 

In  summary,  the  development  of  a  low-cycle  life  prediction  model  that  has  the  capability  to 
account  for  the  effect  of  surface  oxidation  on  life  of  Titanium  matrix  MMC's  is  complete.  The  research 
performed  herein  concluded  that  the  life  of  die  coju^iasite  a|:^ears  tnjhe  eontrolled  by -interface  debonding 
and  subsequent  radial  cracking,  fhese  cracks  appear  to  be  driven  by  a  combination  of  two  factors:  the 
development  of  a  surface  layer  consisting  primarily  of  stoichiometric  TiO?  whiclijBdttees--a''^atational 
eigenstrain;  and  the  embrittlement  of  material  at  the  metal-oxi^'*^  jrtorface.  in  addition,  the  influence  of 
various  heat  treatments  on  the  mechanical  behaviorjif.Sri0  Tf  15-3  MMC  was  also  investigated.  The  study 
revealed  that  the  heat  treatmentseffiseted-ihe  overall  composite  compliance  and  damage  accumulation.  A 
model  was  also  developfid-dttfmg^tiiis  phase  of  the  research,  which  successfully  simulates  the  stress  versijs  '' 
strain  resporLS^^rTQie  composite  accountin^Joriioth  plastic  deformatinn  -md  dam^e.^ 

The  worK  perrormea  under  this  grant  also  included  a  program  to  experimentally  characterize  the 
morphology  of  Ti02,  one  of  the  primary  stoichiometric  oxides  formed  during  oxidation  of  titanium,  in 
order  to  develop  more  accurate  oxide  layer  growth  models.  Ao  part  of  Qii:^  uffou,  IhU  growtn  and  structure 
nf^r  TiOj  mridr  Irnrrr  rrmnniithiii  iininj'ili'  i  iirilt|  ll  [illlily  Htininm  iiii1  Ti 

700^C  fui  tiilie  -lengths  as  short  ao  1  hour  to  a  maximum-of  168  hours.  It  has  been  shown  that  specimen 
geometry  plays  a  significant  role  in  the  rate  of  oxide  growth.  »  Iir~ad4itionj  the  surface  nxjd^tion — 
nf  Ti  1  ^  7  \^Ar  crack  exte^^gi^T^  invpgtifrqtftd-  R^StlltSOf  the 

monotonic  tenfiinn  eypf^rtmPinlQ  nn  prp-rrflrTpH.fhen  oxidized-then  tested  specimens-^dieate  that  the  oxide 
layeTgSLerafiyhas  little  effect  on  the  bulk  behavior  of  the  rnatermt;  pait  oTOus  effort,  nurhefiCkiihodel 
was  also  aeveiopea  which  is  capaole  of  simulating  the  effeef^f  the  oxide  layer  on  the  crack 
resistance  of  metals.  It  is  believed  that  such  a  model  could  be  used  to  facilitate  the  study  of  the  static  and 
fatigue  strength  of  fractured  metal  in  an  oxidizing  environment.  Two  different  methods  of  fracture 
mechanics  modeling  (i.e.  the  isoparametric  and  the  crack  closure  method)  were  evaluated  in  order  to 
simulate  the  mechanical  response  of  a  pre-cracked  pre-oxidized  specimen.  The  current  model  is  capable  of 
describing  the  effect  of  the  stiffer  oxide  and  the  oxide  layer  volumetric  expansion  to  the  energy  release  rate. 

It  has  been  found  that  the  combination  of  these  two  oxide  scale  properties  causes  the  energy  release  rate  ^ 
decrease  when  a  thin  layer  if  oxidf"  is  present  ^nd  to  increag^  for  a  thick  oxide  layer. 

TheTastphagToTthis  researen  effort  was  to  develop  and  numerically  implement  a  mathematical 
model  of  oxidation  for  metals  with  the  capability  of  modeling  complex  oxidation  fronts,  such  as  they  arise 
in  MMC’s.  The  oxidation  of  metals  has  been  modeled  by  modifying  the  Fickian  diffusion  problem  in  order 
to  simulate  the  chemical  reaction  (phase  change)  in  the  metal.  The  current  model  is  c^>able-of  selvmgJjp 
and  2D  oxidation  problems  in  metallic  domains  with  comnleY  eeometry.  Simuiations  of  the  oxide  scale 
gfowtli  from  die  craclr  sui'face~m  a  Ti-15-3  specimen  with  a  semi-mfinite  crack,  a  complex  geometry 
crack-tip,  as  well  as  the  scale  growth  on  cylindrical  and  wedge-like  domain  geometries,  have  been 
successfully  perfbrmedr 
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Status  of  Effort/ Accomplishments/New  Findings: 

The  objectives  of  this  research  effort  were  to:  1)  develop  a  micromechanics  and 
thermodynamics  based  model  for  predicting  macroscopically  averaged  thermomechanical 
constitutive  behavior  of  continuous  fiber  metal  matrix  composites  (MMC’s);  2)  develop  and  test 
MMC’s  experimentally  to  determine  microscopic  variables  that  are  necessary  to  develop  the 
micromechanics  model  mentioned  in  1);  and  3)  carry  out  a  limited  experimental  program  for 
model  verification.  In  the  following  subsections,  key  accomplishments  which  have  lead  towards 
the  fulfillment  of  the  aforementioned  objectives,  are  briefly  summarized. 

Life  Prediction  of  Continuous  Fiber  Metal  Matrix  Composites: 

The  development  of  a  low-cycle  life  prediction  model  that  has  the  capability  to  account 
for  the  effect  of  surface  oxidation  on  life  of  Titanium  matrix  MMC's  is  complete.  We  have 
observed  in  our  own  experiments  (See  Fig.  1)  that  life  is  reduced  by  a  factor  of  five  in  SCS-6/Ti- 
p21S  tested  at  650°C  when  the  composite  was  tested  in  air  rather  than  inert  gas  (Argon).  The  life 
of  the  composite  appears  to  be  controlled  by  interface  debonding  and  subsequent  radial  cracking, 
as  shown  in  Fig.  2.  This  reduction  in  life  appears  to  be  the  result  of  surface  cracks  which  develop 
on  the  specimen  when  tested  in  air,  as  demonstrated  in  the  SEM  photo  shown  in  Fig.  3.  These 
cracks  appear  to  be  driven  by  a  combination  of  two  factors:  the  development  of  a  surface  layer 
consisting  primarily  of  stoichiometric  Ti02  which  induces  a  dilatational  eigenstrain;  and  the 
embrittlement  of  material  at  the  metal-oxide  interface.  The  net  effect  of  these  surface  features  is 
an  overall  reduction  in  the  critical  energy  release  rate  for  the  system.  Using  these  observations 
and  assumptions,  we  have  developed  a  finite  element  model  which  utilizes  a  representative 
volume  element  (RVE)  for  the  laminate  tested  above,  as  shown  in  Fig.  4.  As  described  in  our 
recent  papers,  the  mesh  incorporates  the  surface  oxide  layer,  as  well  as  the  dilatational 
eigenstrain  in  order  to  capture  die  life  reduction  mechanism.  As  shown  in  Fig.  5,  for  the  case  of 
monotonic  loading,  the  model  does  produce  a  reduced  stress-strain  behavior  for  the  composite 
due  to  the  development  of  surface  cracks.  Figures  6  and  7  demonstrate  a  comparison  of  the 
predicted  results  for  cyclically  loaded  unoxidized  and  oxidized  SCS-6/Ti-p21S  [0]4  specimen  at 
650°C,  respectively. 

Influence  of  Heat  Treatment  on  the  Mechanical  Properties  and  Damage  Development  in  a 
SiC/Ti-15-3  MMC: 

The  purpose  of  this  study  was  to  investigate  the  influence  of  possible  heat  treatments  on 
the  mechanical  behavior  of  SiC/  Ti  15-3  MMC,  as  well  as  characterize  the  deformation  and 
damage  mechanisms.  The  effects  of  three  heat  treatment  conditions  on  the  thermomechanical 
response  and  damage  development  in  unidirectional  axial  and  transverse  specimens  were 
investigated.  A  four  ply  uni-directional  SiC/Ti-15-3  MMC  was  tested  in  the  axial  and 
transverse  direction  and  the  damage  evolution  at  both  room  temperature  and  at  427  °C  was 
observed.  The  thermomechanical  study  revealed  that  the  heat  treatments  effected  the  overall 
composite  compliance  and  damage  accumulation. 

In  addition,  averaging  micromechanical  models  were  utilized  to  simulate  the  effect  of 
matrix  damage  and  plasticity  on  the  mechanical  response  of  the  composite.  A  model  based  on 
the  Mori-Tanaka  method  was  developed  and  implemented  which  combines  both  plasticity  and 
damage  in  an  incremental  formulation.  Model  predictions  were  then  compared  to  experimental 
results  through  microstructural  evaluation  of  matrix  crack  densities. 


Mechanical  tests  show  that  the  450°C  heat  treatment  creates  a  microstructure  with  a 
consistently  high  elastic  modulus  and  high  damage  tolerance  for  all  loading  conditions  and 
temperatures.  Microstructural  evaluation  identified  the  primary  damage  modes  for  both  the 
transverse  and  axial  specimens.  The  axial  specimens  showed  evidence  of  cracks  developing 
perpendicular  to  the  loading  direction  starting  from  the  fiber/matrix  interface  (see  Fig.  8).  The 
transverse  specimens  showed  cracks  emanating  from  areas  of  poor  consolidation  resulting  in 
cracks  propagating  in  the  loading  direction  along  grain  boundaries.  The  model  developed 
during  this  research  initiative  successfully  simulates  the  stress  versus  strain  response  of  the 
composite  accounting  for  both  plastic  deformation  and  damage  (see  Fig.  9).  In  addition,  a 
prediction  of  the  crack  density  was  made  as  a  function  of  overall  applied  load.  It  was  shown 
that  the  predicted  crack  density  at  the  final  load  level  over  predicts  the  crack  densities 
measured  from  a  post-test  microstructural  analysis  for  the  transverse  specimens  (see  Fig.  10). 

It  was  believed  that  the  over  prediction  was  most  likely  due  to  pre-existing  matrix/fiber 
damage. 

Experimental  Characterization  of  Oxidation  on  Fracture  Surfaces  and  Crack  Growth  Behavior: 

The  purpose  of  this  research  was  to  experimentally  characterize  the  morphology  of 
Ti02,  one  of  the  primary  stoichiometric  oxides  formed  during  oxidation  of  titanium,  in  order  to 
develop  more  accurate  oxide  layer  growth  models.  In  addition,  this  phase  of  the  research 
included  the  preliminary  design  and  development  of  a  test  methodology  which  will  allow  for 
real  time  evaluation  of  crack  growdi  behavior  in  an  elevated  temperature  (25°to  700°C)  and 
oxidizing  (air)  environment. 

As  part  of  the  effort  to  characterize  the  growth  and  structure  of  the  Ti02  oxide  layer, 
monolithic  samples  of  high  purity  titanium  and  Ti-15-3  were  oxidized  in  air  at  700°C  for  time 
lengths  as  short  as  1  hour  to  a  maximum  of  168  hours.  Various  specimen  geometry’s  were 
oxidized  and  then  analyzed  using  SEM.  Geometry’s  investigated  include:  the  inside  radius  of  a 
small  diameter  hole  («300  pm);  the  outside  radius  of  a  small  diameter  wire  (» 140pm);  and  a 
flat  edge  which  included  a  corner.  In  addition,  in  order  to  gain  insight  on  how  the  oxide  layer 
would  effect  crack  growth  behavior  and  material  fracture  toughness,  a  series  of  experiments 
were  conducted  wherein,  samples  of  high  purity  titanium  were  pre-cracked  to  a  predetermined 
length  using  a  fully  reversed  low  cycle  fatigue  load.  The  crack  in  each  specimen  was  wedged 
open  and  then  either  heat  treated  (in  an  inert  enviromnent)  or  oxidized  at  700°C  for  24  to  48 
hours.  Each  specimen  was  then  subjected  to  a  monotonic  tensile  load  were  visual  and  tabular 
data  were  recorded  for  later  analysis. 

It  has  been  shown  that  specimen  geometry  plays  a  significant  role  in  the  rate  of  oxide 
growth,  as  shown  in  Fig.  11.  However,  additional  insight  was  needed  in  order  to  understand 
the  complex  strain  field  that  appears  to  be  present  in  the  oxide  layer  during  formation,  as  well 
as  the  predominate  diffusion  mechanism.  Figures  12  a)-d)  show  typical  oxide  scale  developed 
during  oxidation  at  700°C.  The  surface  oxidation  characteristics  of  Ti-15-3  MMC  while 
undergoing  crack  extension  is  shown  in  Fig.  12d).  The  photograph  indicates  how  oxide  scale 
from  the  matrix  develops  around  the  exposed  sigma  fibers.  Figure  13  shows  how  an  oxide 
layer  develops  on  the  fracture  surface  during  oxidation.  Results  of  the  monotonic  tension 
experiments  on  pre-cracked-then  oxidized-then  tested  specimens  indicate  that  the  oxide  layer 
generally  has  little  effect  on  the  bulk  behavior  of  the  material  (Fig.  14).  Figure  15  a)-c) 
indicates  that  the  oxide  generally  fractures  at  loads  substantially  below  that  which  are  required 
for  crack  growth  extension  of  the  material  considered  herein. 


Modeling  the  Effects  of  Oxidation  on  the  Crack  Growth  Resistance  of  Metals: 

The  piirpose  of  the  research  was  to  develop  a  numerical  model  capable  of  simulating  the 
effect  of  the  oxide  layer  on  the  crack  growth  resistance  of  metals.  It  was  envisioned  that  such  a 
model  would  then  be  used  to  facilitate  the  study  of  the  static  and  fatigue  strength  of  fractured 
metal  in  an  oxidizing  environment. 

Two  different  methods  of  fracture  mechanics  modeling  (i.e.  the  isoparametric  and  the 
crack  closure  method)  were  evaluated  in  order  to  simulate  the  mechanical  response  of  a  pre¬ 
cracked  pre-oxidized  specimen.  Validation  of  the  methods  is  accomplished  by  comparing 
numerical  results  with  an  analytical  solution  for  a  crack  in  a  homogeneous  material  (see  Fig 
16),  as  well  as  crack-tip  extending  through  a  heterogeneous.  Two  physical  characteristics  of 
the  oxide  scale  (i.e.  the  stiffness  change  and  the  volumetric  expansion)  were  as  input  in  the 
parametric  numerical  studies.  In  particular,  the  effect  of  these  two  parameters  on  the  energy 
release  rate  was  evaluated  for  a  monotonically  loaded,  pre-cracked,  compact  tension  specimen 
at  room  temperature. 

The  current  model  is  capable  of  describing  the  effect  of  the  stiffer  oxide  and  the  oxide 
layer  volumetric  expansion  to  &e  energy  release  rate.  It  has  been  found  that  the  combination 
of  these  two  oxide  scale  properties  causes  the  energy  release  rate  to  decrease  when  a  thin  layer 
if  oxide  is  present  and  to  increase  for  a  thick  oxide  layer,  as  shown  in  Fig.  17.  An  FEM  model 
was  also  used  to  simulate  a  monotonic  tension  test  on  a  CT  style  specimen  that  was 
preoxidized.  Comparison  to  experimental  results  is  shown  in  Fig.  18. 

Modeling  of  Oxidation  Fronts  in  Metals: 

The  purpose  of  this  research  was  to  develop  and  numerically  implement  a  mathematical 
model  of  oxidation  for  metals  with  the  capability  of  modeling  complex  oxidation  fronts,  such 
as  they  arise  in  MMC’s. 

The  oxidation  of  metals  has  been  modeled  by  modifying  the  Fickian  diffusion  problem 
in  order  to  simulate  the  chemical  reaction  (phase  change)  in  the  metal.  The  resulting 
mathematical  model  consists  of  two  parabolic  differential  equations  and  metal-oxide  interface 
conditions.  Two  different  variants  of  a  fixed  grid  finite  element  method  for  numerical 
simulation  of  the  oxidation  process  have  been  used.  The  first  approach  (the  discrete  interface 
method)  taken  was  to  locate  the  oxidation  front  and  split  the  domain  into  metal  and  oxide 
subdomains.  The  second  approach  (the  smearing  front  method)  is  based  on  reformulating  the 
diffusion  equations  in  both  the  oxide  and  metal,  resulting  in  a  single  non-linear  equation  for  the 
whole  domain.  Figures  19  and  20  show  a  comparison  between  numerical  and  analytical  results 
for  oxide  layer  thickness  and  oxygen  concentration  profiles,  respectively. 

The  current  model  is  capable  of  solving  ID  and  2D  oxidation  problems  in  metallic 
domains  with  complex  geometry.  Validated  numerical  models  have  been  used  to  simulate  the 
oxide  scale  growth  from  the  crack  surfaces  in  a  Ti-15-3  specimen  with  a  semi-infinite  crack 
(see  Fig.  21),  a  complex  geometry  crack-tip  (see  Fig  22),  in  addition  to  scale  growth  on 
cylindrical  and  wedge-like  domain  geometries. 
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Cycles  to  Failure 


Figure  2.  Fracture  surface  of  unoxidized  SCS-6/Ti-p21S  [0]4  fatigued  at  482°C  for 
4000  cycles. 


Figure  3.  Surface  crack  propagating  radially  into  the  interior  of  SCS-6/Ti-p21S  [0]4. 
Confirmation  of  embrittled  layer. 
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Figure  5.  Comparison  of  elastic  and  viscoplastic  analysis  for  monotonic  loading  of 
SCS-6/Ti-p21S  [0]4  at  650°C. 
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a)  As-Received 


b)  450°C  for  24  hours 


c)  700°  C  for  24  hours 


Figure  8.  Etched  microstructure  in  transverse  specimens  of  SiC/Ti-15-3  for  each  heat 
treatment  showing  dominant  crack  growA  along  grain  boundary. 
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Figure  9.  Stress  vs.  strain  for  transverse  specimen  of  SiC/Ti-15-3  comparing 
simulation  utilizing  the  crack  density  prediction  and  experimental  data. 
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Figure  10.  Prediction  of  crack  density  to  applied  stress  for  SiC/Ti-15-3 
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Figure  11.  Oxide  layer  thickness  as  a  function  of  time  for  various 
geometries  of  Ti  at  700°C. 


a)  Wire  Ti  specimen  oxidized  for 
168  hours. 


c)  Hole  in  Ti-15-3  specimen 
oxidized  for  48  hours. 


b)  Flat  corner  of  Ti-15-3  specimen 
oxidized  for  48  hours. 


d)  4  ply  SiC/Ti-15-3  CT 
specimen  oxidized  for  10 
horns  at  700°C. 


Figure  12.  Various  oxidized  geometries  at  700°C. 


Figure  13 .  Macro  crack  in  Ti  that  has  been  oxidized  for  24  hours. 


Oxidized  Ti  at  700C 
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Figure  14.  Crack  opening  displacement  as  a  function  of  load  for  an 
oxidized  Ti  specimen.  The  points  a,b,  and  c  correspond 
to  loading  conditions  of  30,  50,  and  70  Lb.’s 
respectively  as  shown  in  Fig.  15. 
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Figure  15.  a)Fatigue  crack  in  Ti  subjected  to  a  30  Lb.  load  (point  a  in  Fig.  14) 
Exposed  surface  (at  left)  and  oxidized  layer  showing  (at  right). 
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Figure  15.  b)  Fatigue  crack  m  Ti  subjected  to  a  50  Lb.  load  (point  b  m  Fig.  14) 
Exposed  surface  (at  left)  and  oxidized  layer  showing  (at  right). 
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Figure  15.  c)  Fatigue  crack  in  Ti  subjected  to  a  70  Lb.  load  (point  c  in  Fig.  14).  Exposed 
surface  (at  left)  and  oxidized  layer  showing  (at  right). 


Figure  16.  Simplified  oxide  geometry  for  the  parametric  studies 
assumes  the  crack  to  be  planer  and  the  oxide  scale 
to  be  surrounded  by  unoxidized  metal. 
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Figure  17.  Numerical  results  showing  the  energy  release  rate  calculated 

inside  the  oxide  layer  by  the  isoparametric  method  and  the  crack 
closure  method. 
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Figure  18.  FEM  model  used  to  simulate  monotonic  test  on  an  oxidized  Ti 
CT  specimen  (above)  and  COD  vs.  load  comparison  (below). 
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Figure  20.  Comparison  of  oxygen  concentration  profiles  along  X2  perpendicular  to  the 
crack  surface. 


Figure  21.  Simulation  of  the  location  of  the  oxidation  front  in  an  oxidized  Ti-15-3 
specimen  at  700°C  after  0.25  hours,  0.5  hours,  and  1  hour. 


Figure  22.  Simulation  of  crack  tip  in  oxidized  Ti.  Actual  crack 
tip  shown  above  with  simulation  below  it. 


